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摘要 : 为 研究 VLBI 对 探测 器 的 数据 处 理 方法 并 评估 其 处 理 能 力 ， 现 有 做 法 需 基 于 测 站 
VLBI 终端 接收 的 探测 器 信号 进行 。 采 用 信号 仿真 方法 可 以 根据 设计 轨道 与 信 标 特点 ， 利 用 
计算 机 生成 需要 的 探测 器 VLBI 信号 ， 相 比试 验 观测 具有 独特 的 优势 。 观 测试 验 时 ， 测 站 终 
端 接 收 探测 器 下 行 射频 信号 ， 通 过 混 频 变换 为 中 频 信 号 ， 再 进行 数字 化 采集 数据 。 仿 真 信号 
时 ， 为 减少 计算 规模 ， 直 接 构造 数字 中 频 信 号 ， 二 次 采样 后 提取 通道 信号 ， 获 得 数字 终端 的 
VLBI 仿真 信号 。 由 于 探测 器 相对 VLBI 测 站 运动 ， 测 站 接收 的 探测 器 下 行 信号 反映 目标 的 视 
; 向 速度 变化 。 根 据 探 测 器 信号 时 延 模型 ， 研 究 了 探测 器 点 频 信 号 和 有 限 带 宽 信 号 多 普 勒 效 应 
了 的 模拟 方法 。 通 过 仿真 信号 与 VLBI 终端 接收 信号 的 对 比 以 及 对 VLBI 相关 处 理 结果 的 分 析 ， 


验证 了 探测 器 VLBI 信号 仿真 方法 的 可 行 性 ， 为 后 续 仿真 技术 在 月 球 与 火星 探测 器 中 的 VLBI 
测定 轨 应 用 黄 定 了 技术 基础 。 
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其 长 基线 干涉 测量 "1( Very Long Baseline Interferometry, VLBI) 是 20 世纪 60 年 代 后 期 出 现 的 一 种 

射电 天 文 技 术 。 它 具有 空间 分 辨 率 高 、 测 量 精度 高 的 特点 ， 现 已 广泛 应 用 于 天 体 物理 、 天 体 测量 、 天 
文 地 球 动力 学 和 深 空 探测 领域 。 我 国 于 2007 年 首次 在 探 月 工程 中 使 用 VLBI 技术 ,建立 了 中 国 
VLBI 网 ( Chinese VLBI Network, CVN)， 包 括 上 海 余 山 站 、 北 京 密云 站 、 云 南 昆 明 站 、 乌 鲁 木 齐 南山 
站 和 上 海 VLBI 数据 处 理 中 心 。 中 国 VLBI 网 已 参加 了 探 月 工程 娣 娥 系列 探测 器 的 多 次 任务 ， 对 关键 
弧 段 的 快速 精密 测定 轨 发 挥 了 不 可 替代 的 重要 作用 1。 
i VLBI 测 站 由 天 线 、 接 收 机 、 数 据 采 集 终端 、 时 间 频 率 系统 等 组 成 ， 工 作 时 ， 各 测 站 接收 机 天 线 
© 汇集 的 射频 信号 放大 后 ， 经 独立 本 振 进行 混 频 得 到 中 频 信 号 ， 然 后 传送 至 数字 化 终端 ， 高 速 模 数 采 样 
后 作 数字 基带 转换 ， 得 到 所 需 频 率 通 道 的 基 频 信号 ， 最 后 经 过 编码 ， 记 录 在 磁盘 或 者 直接 通过 网 络 传 
送 到 VLBI rat», 

为 研究 VLBI 中 心 对 不 同 探测 需 的 测量 方法 ， 评 佑 分析 中 心 各 功能 模块 的 能 力 及 整个 系统 的 正确 
性 、 可 靠 性 ， 目 前 常规 做 法 是 在 每 次 飞行 任务 前 开展 大 量 实 际 目标 联 测 。 对 单 目 标 探测 器 的 实际 观测 
演练 可 评估 测 站 与 系统 的 可 靠 性 。 但 对 于 那些 全 新 的 深 空 探测 任务 ， 在 正式 任务 前 有 时 难以 找到 合适 
的 探测 器 作为 试验 对 象 。 比 如 在 后 续 探 月 工程 中 ， 将 要 开展 针对 月 球 轨道 的 交会 对 接 ， 届 时 需要 对 月 
球 距 离 的 两 个 动态 目标 进行 同 波束 测量 。 显 然 ， 目 前 没有 类 似 的 实际 探测 器 可 作为 动态 双 目 标的 观测 
对 象 。 若 采用 探测 器 VLBI 信号 仿真 方法 ， 则 可 根据 探测 器 的 数量 、 轨 道 、 信 标 等 特点 ， 模 拟 生 成 
VLBI 终端 接收 的 不 同 探测 器 信号 ， 提 供 月 球 轨道 双 目 标 同 波束 VLBI 测量 数据 ， 用 于 研究 和 改进 
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VLBI 处 理 模块 。 通 过 仿真 与 实际 演练 两 种 方法 结合 ， 可 较 好 地 评 佑 VLBI 中 心 与 整个 系统 对 双 ( 多 ) 
目标 探测 器 的 测量 能 力 。 因 此 ， 研 究 VLBI 信号 仿真 技术 有 重要 意义 与 迫切 工程 需求 。 

在 探 月 工程 中 ， 探 测 需 下 行 VLBI 射频 信号 为 S/X 频段 ， 采 用 调 相 方式 ， 调 制 信号 在 形式 上 主要 
包括 点 频 信号 和 有 限 带宽 信号 两 类 。 在 仿真 信号 时 ， 没 有 直接 仿真 射频 信号 ， 而 是 先 将 信号 从 射频 变 
换 至 数字 中 频 ， 然 后 再 对 中 频 信号 采样 ， 提 取 通 道 信 号 ， 最 终生 成 用 于 验证 VLBI 数据 处 理 系统 的 仿真 
数据 ， 这 样 既 完 整地 保留 了 调制 信号 频谱 的 主要 特性 ， 又 避免 了 直接 处 理 需 极 高 采样 率 的 射频 信号 。 


1 探测 带 信 号 结构 


在 探 月 工程 中 ， 探 测 器 VLBI 下 行 射 频 信号 采用 调 相 方式 ,包括 遥测 信号 、 单 向 差分 测 距 
(Differential One-way Ranging, DOR ) 信 号 、 测 距 信 号 、 残 留 遥 控 信 号 (图 1)。 其 中 ,下 行 载波 、 测 距 
证 号 ( 测 距 主 音 和 测 距 次 音 ) 和 单 向 差分 测 距 信 号 为 点 频 信 号 ; 遥测 信号 和 遥控 信号 为 经 过 相位 调制 
的 有 限 带宽 信号 。 这 些 信号 统一 调制 在 下 行 信号 的 主 载波 上 ， 数 学 模型 ”为 


s(t) = /2P, sin[ 2mf.t + 0,d (t) * 0,d,(t) + d,sin(2mfit) + d,sin(2mf,t) Ie (1) 
~ 其 中 ，P; 为 总 功率 ; f 为 载波 频率 ; d (让 为 遥测 信号 ; 0, 为 遥测 信号 的 调制 指数 ，d (2) 为 遥控 信 
一 号 ; 0, 为 遥控 信号 的 调制 指数 ; /, 为 测 距 信号 ;中 , 为 测 距 信号 的 调制 指数 f 为 单 向 差分 测 距 信号 ; 
中 为 单 向 差分 测 距 信 号 的 调制 指数 。 
= 下 行 载波 
残留 遥控 测 距 次 音 遥测 测 距 次 音 MEES DORY 1 DOR Z2 


8 kHz 16~20 kHz 65.536 kHz 100 kHz 500 kHz 3.8 MHz 19.2 MHz 
fI f 


图 1 X 波段 下 行 信号 上 边 带 频谱 示意 图 
Fig.1 The spectrum of X band downlink signal 


2 仿真 原理 与 实现 方法 


2.1 探测 需 信 号 时 延 模型 构造 

VLBI 测 站 与 探测 需 之 间 存 在 相对 运动 ， 因 此 在 仿真 时 要 考虑 相对 运动 效应 引起 的 测 站 接收 探测 
右 信 号 的 时 间 延 迟 (时 延 ) 及 其 变化 沾 。 时 延 包含 几何 时 延 、 引 力 时 延 、 天 线 时 延 修 正 、 源 结构 时 
延 修 正 、 中 性 大 气 与 电离 层 时 延 修正 等 ， 其 中 ， 几 何 时 延 是 理论 时 延 的 主要 部 分 ， 随 目标 的 运动 
而 变化 。 

设 地 心 0、 测 站 4 及 探测 器 B( 图 2)，g(i) 表 示 ;i 时 刻 参考 点 (地 心 ) 接 收 到 信号 波 前 的 光 行 时 ; 
gi(i) 表 示 1 时 刻 参 考点 (地 心 ) 信 号 波 前 超前 到 达 测 站 的 光 行 时 ; es Ct) 表示 i; 时刻 参考 点 ( 测 站 ) 接 收 
到 信号 波 前 的 光 行 时 。 若 妃 时 刻 从 探测 器 她 发 出 信和 号， 同一 波 前 到 达 地 心 0、 测 站 A 的 时 刻 分 别 为 
tos ti, MARIER, goo), gio) A g(t) DIA 
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go(fo) = to — t5, 


gi(49) 7 tj 7 tos (2) 
g(h)7t -ts 
显然 ， 
82(t1) = go(to) + gi(to). (3) 
在 实际 应 用 中 ， 每 个 观测 弧 段 的 时 延 模型 都 以 五 次 时 
延 多 项 式 表示 : 


g(t) 2b, + but t b, +b + bt + bu, (4) 


FER, by, bi, .… ,bs 为 时 延 模型 多 项 式 系数 。 

根据 实际 时 延 模型 g(1) ，&i(2) ， 用 最 小 二 乘法 作 五 
次 多 项 式 拟 合 ， 得 到 时 延 模型 g,(1) 。 
2.2 时 延 模型 线性 近似 

为 减少 计算 量 ， 在 短 时 间 内 假设 时 延 模型 随时 间 线 性 


变化 ， 即 

g'a(t) = b, + bit, (5) 图 2 VLBI 时 延 示意 图 
其 中 ， RA bo, b, 在 短 时 间 内 为 恒 常数 ， 忽 略 时 延 模型 g， xs 
(1) 高 次 项 。 


在 短 时 间 内 用 (5) 式 表示 线性 化 时 延 模型 存在 一 定 误差 ， 误 差 大 小 取决 于 用 线性 近似 的 时 间 长 度 
和 五 次 时 延 多 项 式 中 高 次 项 系数 。 相 关 人 处 理 机 '" 使 用 皮 秒 量 级 的 时 延 模型 精度 足以 满足 应 用 需求 ， 
故 误 差 值 在 0.01 ps(1e-14s) 以 内 的 模型 误差 即 可 视 为 满足 仿真 精度 要 求 。 在 较 短 的 观测 时 间 内 ， 时 
延 模 型 g (t) 的 二 次 项 是 线性 近似 的 主要 误差 项 ， 其 系数 b, 一 般 为 e-11 量 级 ， 要 想 使 得 该 项 引起 的 
时 延误 差 小 于 0. 01 ps， 选 定 的 线性 化 时 间 长 度 应 小 于 0.01s。 其 他 高 次 项 的 系数 5,, b, bs, JRE 
e-14 量 级 ， 由 此 引起 的 误差 可 忽略 。 
2.3 构造 中 频 信和 号 

探测 器 下 行 载波 信号 的 类 型 主要 包括 点 频 信 号 和 有 限 带宽 信号 。 其 中 点 频 信和 号 包括 : 载波 、 测 距 
信号 和 单 向 差分 测 距 信号 。 有 限 带 宽 信 号 包括 遥测 信号 和 遥控 信号 ， 是 以 正弦 波 为 载波 的 调 相 信和 号 ， 
调制 方式 为 PCM/PSK/PM， 一 般 可 表示 为 sin(2mwfi+0)， 用 相位 9 的 变化 传递 数字 信息 。 对 于 二 进 制 
相 移 键 控 信 号 ，6 以 0/m 交替 变化 表示 二 进 制 位 0/1, 1 对 应 相位 0，-1 对 应 相位 m. DRE DU RUE TE 
ESLI RIRN Asin(2mfi) , A 随 相位 0 的 变化 而 变化 。 

为 完整 地 保留 调制 信号 及 其 边 带 频谱 ， 选 择 中 频 信和 号 约 为 60 MHz。 有 限 带 宽 信 号 的 码 速率 最 大 
为 2 048 bits， 码 元 变化 周期 是 采样 周期 的 两 万 九 千 多 倍 。 在 一 个 码 元 周期 内 ， 有 限 带宽 信号 保持 不 
变 ， 对 其 采样 可 仿真 遥测 和 遥控 有 限 带宽 信号 。 因 此 ， 需 选择 合适 的 时 间 尺 度 对 点 频 信 号 和 有 限 带宽 
信和 号 采样 ， 仿 真 中 频 信和 号 。 

为 简约 起 见 ， 根 据 (1) 式 探测 器 转发 的 下 行 信号 的 一 般 形式 ， 分 析 其 中 典型 的 点 频 信号 、 有 限 带 
宽 信 号 多 普 勒 效应 影响 。 信 号 表达 式 如 下 .: 


f(t) = sin[ 2uf.t + OAsin(27ft) + 9,sin(27hi) ], (6) 


其 中 , f. ,分别 为 载波 和 点 频 单 向 差分 测 距 信号 ; f 为 有 限 带宽 遥测 信号 ; 4 为 遥测 信号 的 幅 值 ; 
01 ,中 , 为 调制 系数 。 


auos c 时 刻 接收 的 信号 (i) 是 探测 器 在 i-g,(1) 时刻 发 出 的 ， 短 时 间 内 假设 时 延 模型 随时 间 线 性 


变化 , 用 g CO TRE g,(1) ，f(1) 与 探测 器 转发 的 f(1) 关 系 是 
JO) =fli - 2,0], (7) 
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将 (5) 式 和 (6) 式 代入 (7) 式 ， 得 到 : 
FG) 2 sin(2mf.[ (1 — b,)t - b] * 6,Asin2mf, [ (1 — b,)t - b] 
+ ó,sin2mf;[ (1 -b,)t -b,]}. 

对 比 (6) 式 ， 测 站 接收 信号 (4) 的 载波 /. 和 信号 户 , fs 发 生 频 率 偏 移 ， 体 现 了 信号 传输 中 由 于 相对 运 
动产 生 的 多 普 勒 效应 。 频 率 偏 移 量 分 别 为 bf Fl b Cf. bU E Ee 

根据 (8) 式 ， 探 测 器 下 行 信 号 在 频 域 上 表现 为 关于 载波 左右 对 称 ， 其 频谱 主要 分 布 在 [ (1-5) 
eau). Cleh et fr) ls HP fr N max, fa), Klo 

构造 中 频 信 号 f(1) 如 下 : 


^ 


f) = sinQm[f.(1 7 6) -falt -fbo} + 0,Asin2mf, [ (1 - 5,)t - bo] 
十 $,sin2mf;[ (1 = b)i = b, ] , 


(8) 


(9) 


对 比 (9) 式 与 (8) 式 ， 只 是 载波 频率 发 生变 化 ， 其 余 信 号 不 变 。 即 通过 降低 载波 频率 ， 将 信号 从 
射频 搬移 到 中 频 ， 且 信号 方 , f 前 后 保持 不 变 。 

综 上 ，VLBI 测 站 接收 信号 是 添加 时 延 的 射频 信号 ， 由 (1) 式 添加 时 延 并 降低 载波 频率 ， 得 到 包 
含 多 普 蔓 效应 的 中 频 信号 ， 形 式 如 下 : 


s(t) = /2P, sin(2m[f.(1 - b1) -falt - f by) 
+ 6, Asin2uf, [ (1 - b,)t - b,] * 6,A'sin2mf [ (1 — b,)t — b] (10) 
+ $,sin2nmf[(1— b,)t - bi] + ,sin2mf;[ (1 - b,))t - bo], 
HP, fA 为 遥控 信号 ; 4 为 遥控 信号 幅 值 。 
2.4 VLBI 终端 信号 仿真 
图 3 显示 了 仿真 VLBI 终端 接收 信和 号 实现 流程 。 先 仿真 带 有 多 普 
勒 效 应 的 探测 器 VLBI 下 行 中 频 信和 号， 再 滤波 得 到 某 个 通道 的 频率 区 
间 ， 通 过 希 尔 伯 特 变换 保留 正 频率 部 分 ， 再 通过 数字 变频 移 频 至 基 
带 ， 经 希 尔 伯 特 逆 变 换取 实 信号 ， 最 后 通过 重 采样 降低 采样 率 得 到 通 
道 实 信 和 号。 通道 宽度 为 2 MHz 时 ,采样 率 取 4 MHz, 


3 数值 试验 


3.1 仿真 探测 器 下 行 中 频 信号 

根据 上 述 方法 采用 某 月 球 探测 器 信 标 实际 参数 仿真 了 探测 器 
VLBI 下 行 中 频 信和 号， 其 中 包含 测 距 信 号 (500 kHz EFM). EM 
信和 号 ( 副 载 波 65. 536 kHz 且 码 速率 256 bit/s) 和 两 对 单 向 差分 测 距 信 
号 。 高 频 信 号 起 始 频 点 了 为 8460 MHz， 降 低 载波 频率 后 生成 探测 器 
下 行 中 频 信号 ， 采 样 率 117 MHz， 快 速 侍 里 叶 变 换 点 165 000。 

图 4 显示 带宽 58.5 MHz 的 探测 器 下 行 中 频 信号 频谱 。 图 中 ， 箭 图 3 信号 仿真 流程 图 
头 3 为 载波 信号 ，2 和 4 是 距 载波 约 3. 85 MHz 的 第 1 对 单 向 差分 测 距 信 Fig. 3 Steps of signal simulation 
号 (DOR1); 1 和 5 是 距 载 波 约 19. 27 MHz 的 第 2 对 单 向 差分 测 距 信 号 
(DOR2) ， 载 波 两 侧 500 kHz 处 是 测 距 主 音信 和 号。 另外， 有 限 带宽 的 遥测 信号 (带宽 约 0. 5 kHz) 靠近 载波 。 

测 距 信 号 、 琐 测 信 号 、 单 向 差分 测 距 信 号 关于 载波 对 称 ， 调 制 系数 决定 各 信号 的 功率 分 配 ， 其 中 
载波 、 测 距 信 号 和 单 向 差分 测 距 信 号 较为 明显 。 由 于 快速 侍 里 叶 变 换 的 栅栏 效应 使 图 中 载波 两 侧 对 应 
频 点 幅 值 产生 细微 差别 。 其 他 信号 谐 波 项 能 量 较 小 ， 幅 度 基本 在 -180 dB 以 下 。 
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图 4 探测 器 下 行 中 频 信号 频谱 图 


Fig.4 The spectrum of downlink probe quasi-intermediate-frequency signal 


3.2 ”可行 性 验证 
通过 仿真 信号 与 VLBI 终端 接收 信号 的 对 比 以 及 对 VLBI 相关 处 理 结果 的 分 析 ， 验 证 了 探测 器 
VLBI 信号 仿真 方法 的 可 行 性 。 

图 5(a) 为 实际 信和 号 单 通道 自 功率 谱 ， 从 中 可 观察 到 0. 96 MHz 处 的 载波 信号 ; 距离 载波 500 kHz, 
65. 536 kHz 的 测 距 主音 信号 和 遥测 信号 ; 其 他 为 信号 谐 波 项 。 

图 5(b) 为 仿真 VLBI 终端 接收 信号 自 功 率 谱 (采样 率 4 MHz, 快速 侍 里 叶 变 换 点 262144， 采 样 
1 s) ， 包 含 载波 信号 、 测 距 信 号 (500 kHz 主音 频率 ) 、 遥 测 信和 号 ( 副 载波 65.536 kHz 且 码 速率 1 024 
bit/s) 、 高 斯 白 噪声 信和 号。 频率 0. 968 MHz 处 为 载波 信号 ， 测 距 主 音信 号 和 有 限 带 宽 的 遥测 信号 分 别 
距离 载波 500 kHz, 65.536 kHz， 关 于 载波 对 称 。 距 离 测 距 主 音 65. 536 kHz, 131.072 kHz 且 关 于 测 距 
主音 对 称 的 为 遥测 信号 谐 波 项 。 
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(a) (b) 
图 5 单 通道 自 功率 谱 图 (a) 实际 信号 ; (b) 仿 真 信号 
Fig. 5 The auto-power spectrum of channel signal 


(a) The auto-power spectrum of real signal; (b) The auto-power spectrum of simulation signal 


仿真 信号 与 对 应 的 实际 信号 对 比 表 明 ， 载 波 信号 、 测 距 信 号 、 遥 测 信号 的 频率 和 功率 分 配 与 实际 
信号 一 致 ， 噪 声 信号 的 功率 和 分 布 也 与 实际 信号 相同 ， 验 证 了 仿真 方法 的 可 行 性 。 
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探测 器 VLBI 下 行 中 频 信 号 分 别 添加 不 同时 延 模型 ， 仿 真 两 台 站 终端 接收 信号 (采样 率 4MHz， 快 
速 傅 里 叶 变换 点 512， 采 样 1 s) 。 经 过 VLBI 相关 处 理 ， 可 得 互 功率 谱 帆 度 和 相位 图 (图 6)。 图 6(a) 的 
幅度 谱 显 示 ， 载 波 信号 和 遥测 信号 在 [0.8 MHz, 1.2 MHz] 频段， 遥测 信号 的 谐 波 项 和 测 距 信和 号 在 
[0.4 MHz, 0.6MHz]，[1.4MHz,，1.6MHz] 频 段 。 图 6(b) 的 相位 谱 显 示 ， 在 载波 信 导 、 遥 测 信和 号、 
测 距 信 号 分 布 的 主要 频段 ， 相 位 谱 有 明显 相关 条 纹 ， 其 它 频段 的 点 由 于 信号 较 弱 ， 相 位 分 布 弥散 ， 条 
纹 较 弱 ， 验 证 了 不 同 台 站 接收 的 探测 器 VLBI 下 行 仿真 信号 的 相关 性 。 
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图 6 仿真 两 站 信号 互 功率 谱 图 (a) 幅 度 谱 ; (b) 相位 谱 与 条 级 


Fig. 6 The cross-power spectrum of two station simulation signal. (a) The amplitude spectrum; (b) The phase spectrum and strip 
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本 文 在 研究 探测 器 信号 结构 的 基础 上 ， 分 析 了 包含 探测 器 下 行 点 频 信 号 和 有 限 带宽 信号 的 VLBI 
终端 信号 仿真 方法 ， 构 造 了 中 频 信号 仿真 VLBI 数字 终端 接收 信号 。 采 用 最 小 二 乘法 可 构造 出 满足 精 
度 的 时 延 模 型 ， 用 于 模拟 多 普 勒 效应 。 通 过 仿真 数字 终端 信号 与 实际 信号 的 比较 以 及 对 不 同人 台 站 的 仿 
真 信号 作 相关 处 理 ， 验 证 了 探测 器 VLBI 信号 仿真 方法 的 可 行 性 。 后 续 建 立 仿真 各 模块 误差 分 析 模 
型 ， 完 善 仿 真 系统 ， 为 月 球 探测 器 VLBI 信号 仿真 的 应 用 ， 疙 至 火星 探测 器 VLBI 信号 仿真 的 研究 打 
下 良好 的 基础 。 
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Abstract; VLBI ( Very Long Baseline Interferometry ) technique is an important method for orbit 
measurement of deep space probe. In order to research the data processing methods of the VLBI center and 
evaluate the measurement capability of the VLBI system, at present the routine method is observing actual 
deep space probe by carrying out large observation experiment to get probe signal that is received by VLBI 
terminals. Compared with actual observation, simulating probe signal based on the characteristics of the track 
and beacon is a better choice. The actual downlink signal is radio-frequency signal, which will be converted to 
intermediate-frequency signal, then collected by receiver. In simulation, we construct the intermediate- 
frequency signal directly to reduce the amount of calculation, and generate the digital terminal signal via 
resampling. Because of the relative motion between the probe and the VLBI station, the frequency of the signal 
shifts accordingly, which is called Doppler effect. To simulate the actual probe signal, we build the delay 
model, and investigate the simulation method about the Doppler effect of the typical point frequency signal and 
finite bandwidth signal. By comparing the simulated signal with the actual signal received by VLBI terminal, 
and analyzing the VLBI correlation result, the feasibility of the simulation method is verified, which paves the 
road for the following VLBI application in the moon and Mars. 

Key words: Signal simulation; Deep space exploration; VLBI; Beacon 


